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Abstract 

We consider a two-user multiple-input multiple-output (MIMO) interference channel (IC), where a 
single data stream is transmitted and each receiver applies the minimum mean square error (MMSE) 
filter. In this paper, we study an open topic on the Pareto boundary of the rate region. The Pareto 
boundary is divided by two turning points into the weak Pareto boundary (including the horizontal part 
and vertical part) and the strict Pareto boundary (including the upper-right part and turning points). The 
weak Pareto boundary and turning points can be computed exactly. For the strict Pareto boundary, we 
propose a computationally efficient method called iterative alternating algorithm (lAA) for maximizing 
the rate of one user while the rate of the other user is fixed. To deal with the difficult coupling of the two 
transmit beamformers in this optimization problem, we convert it into two single-beamformer optimization 
problems. Then, by certain equivalent transformations, each problem becomes a quadratically constrained 
quadratic programming (QCQP) problem, which can be optimally solved by the semidefinite relaxation 
(SDR) technique in combination with a matrix rank-one decomposition. Furthermore, convergence of the 
proposed lAA is guaranteed. Numerical simulations show that the lAA provides a better lower bound 
on the strict Pareto boundary compared with the existing methods. 



I. Introduction 

In cellular wireless systems, multiple sectors of different cells share the same time-frequency resource 
for transmission in order to increase the spectral efficiency and occupancy level. Obviously, inter- 
cell or inter-sector interference creates the interference channel (IC). For medium and high SNRs and 
transmitters/receivers operating close to the cell edge, the IC is interference limited. In this paper, we 
apply the achievable rate of each Unk as the performance measure and assume single-user decoding (SUD) 
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at the receivers. The outermost boundary of the achievable rate region is called the Pareto boundary. For 
Pareto-optimal points, it is impossible to increase the rate of one user without decreasing the rate of the 
other user simultaneously. Thus, it is desired to design resource allocation schemes to achieve Pareto- 
optimal points. It has been shown that the competitive operation of the links leads to Nash equilibrium, 
which is usually far away from the Pareto boundary |1|. The coordination between the transmitters is 
required in order to improve their joint outcome and guarantee efficient system operation. 

How to design cooperation schemes for the IC to achieve the Pareto boundary has attracted intensive 
research for several decades. Recently, some important results on characterization or computation of 
the Pareto boundary for the multiple-input single-output (MISO) IC have been obtained, and also the 
multiple-input multiple-output (MIMO) IC has been studied under some special conditions. Here, we 
give a brief, comprehensive, yet incomplete review of cooperation schemes for the MISO/MIMO IC. 

A. MISO IC 

For a two-user MISO IC, all possible Pareto-optimal beamformers can be characterized as linear 
combinations of the zero-forcing (ZF) and maximum-ratio transmission (MRT) beamformers |2|. This 
characterization is later used to derive a closed-form one-dimensional expression for the achievable rate 
region in 13). A general framework for parameterizing Pareto-optimal transmit strategies was proposed 
in which is applicable when the utility functions of the systems are monotonic in the received power 
gains. From the perspective of optimization, computing the Pareto-optimal points may boil down to solving 
a sequence of feasibility problems after certain transformations, which can be converted into second- 
order cone programming (SOCP) problems. This idea motivates an efficient algorithm for maximizing 
the SINR of one link for a given fixed SINR of the other link in a two-user MISO IC fSl. A decentralized 
algorithm based on the iterative updates of certain interference temperature constraints across different 
pairs of transmitters was proposed in [6||, and two distributed optimal algorithms based on alternating 
projections for the MISO IC were developed in Q. 

B. MIMO IC 

However, it is not straightforward to extend the methods for the MISO IC to the MIMO IC, because 
the rate (or SINR) of the MIMO IC depends on both transmit and receive strategies. The hard coupling of 
transceiver strategies makes it extremely difficult to exactly achieve the whole Pareto boundary. Therefore, 
to the best of our knowledge, current studies on the MIMO IC have focused on efficiently finding high 
quality suboptimal operating points. 
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A commonly adopted method to obtain the Pareto boundary for the MIMO IC is via solving a sequence 
of weighted sum-rate maximization problems. Generally, this problem is non-convex and NP-hard, even in 
the single antenna case 18]. Therefore there is no guarantee to obtain a global optimal solution. A linear 
transceiver design algorithm for weighted sum-rate maximization for the MIMO interfering broadcast 
channel (IBC) was recently proposed in |9|, where the sum-rate maximization problem is converted into 
an equivalent sum-MSE cost minimization problem. The iterative algorithm is guaranteed to converge to 
a stationary point. In addition, there also some weighted sum-rate maximization problem with additional 
constraints. For instance, the interference-pricing concept ifTOl is applied to maximize the sum rate of 
the two-user MIMO IC in 1 1 1 1 ; Combining the parameterization lH and distributed interference pricing, 
a beamforming technique was proposed based on balancing the egoistic and the altruistic behavior to 
maximize the sum rate for the single beam MIMO IC in |[T2l : In |[T3l . the authors also studied cooperative 
algorithm to maximize sum-rate for the multi-user MIMO IC based on concept of interference alignment 
|[T4l . However, it is known that |[T4l well studied the degree of freedom of the MIMO IC rather than 
capacity region. 

The weighted sum-rate maximization achieves a sequence of maximum weighted sum-rate operating 
points with different weights on the convex hull of the achievable rate region. However, this approach 
cannot guarantee the finding of all the Pareto boundary points when the rate region is non-convex because 
the tangent lines cannot be achieved at the non-convex part of boundary (S]. Furthermore, when there 
exists a rate constraint, it is not clear how to obtain the corresponding weights to achieve a rate tuple by 
the weighted sum-rate maximization. 

A few works study the achievable rate region in the MIMO IC from other perspectives. For example, 
two interference aware-coordinated beamforming (lA-CBF) algorithms for the two-user MIMO IC were 
proposed in |[T5l . However, the MMSE lA-CBF achieves a lower bound of the sum rate, and the ZF 
lA-CBF achieves some operating points in the rate region by ZF strategies. Recently, a parametrization 
to characterize the Pareto boundary of the multi-cell MIMO performance region was proposed in llT6l . 
However, under an assumption that each receiver (so called simple receiver) has only a single effective 
antenna in the optimization, it clearly leads to performance-loss on the achievable rate region comparing 
with joint transceiver optimization. 

C. Outline and Contributions 

In this paper, we consider a two-cell environment, where each cell has a base station (BS) with 
multiple antennas and a mobile station (MS) with multiple antennas. This scenario can be considered as 
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a special case of the interfering broadcast channel (I-BC) for the downlink or the interfering multi-access 
channel (I-MAC) for the uplink, which can be modeled as two MIMO interfering links, TXi i— )• RXi and 
TX2 I—)- RX2, where TXj, i = 1,2 and RXj are i-th transmitter and receiver, respectively. A single data 
stream is transmitted. We assume that the RXs treat interference as noise and that the TXs have perfect 
local channel state information (CSI) and cooperate with each other. Each receiver applies a MMSE filter. 
For this two-user single-beam MIMO IC, how to achieve the Pareto boundary remains an open problem. 

When one link has a rate/SINR constraint (e.g. rate requirement of primary link in underlay cognitive 
radio environments or for private messages), there is a need to find a way to directly compute a Pareto- 
optimal point corresponding to the rate constraint. This need motivates us to propose an optimization 
problem to maximize the rate of one link for a given rate of the other link. We develop an efficient 
algorithm to compute the Pareto boundary for this setting. 

The main contributions to this open problem in the paper are as follows: 

1) We propose another form of the SINR expression based on the Hermitian angle (Proposition 1 in 
Section II-B), which is helpful to understand and analyze the coupling of the transmit beamformers. 

2) We prove that the strict Pareto boundary can be only achieved when both transmitters spend full 
power (Proposition 2 in Section III-A), which is a necessary condition for the strict Pareto-optimal 
beamformers. 

3) Based on Proposition 1-2, we compute exactly some key points: all the weak Pareto-optimal points, 
the turning points (exactly on the strict Pareto boundary) and some of the ZF points (Section III-B). 

4) We propose the lAA to compute the strict Pareto boundary. The problem of maximizing one rate with 
a constraint of the other rate is a non-convex optimization problem. Firstly, this problem is separated 
into two single-beamformer optimization problems by fixing each transmit beamformer, respectively. 
Then, by certain equivalent transformations, each single-beamformer optimization problem can be 
recast as a QCQP problem, which can be optimally solved by the SDR technique ifTTl . ifTSl in 
combination with a matrix rank-one decomposition [19]. The convergence of the lAA is guaranteed 
(Section IV-B). 

D. Notation 

= and =^ stands for "is defined as" and "perform eigen-decomposition as", respectively. C*^^"* denotes 
the nxm complex matrix space. (•)*, (•)^, (•)^ and (•)^ denote complex conjugate, transpose, Hermitian 
transpose and Moore-Penrose pseudo inverse, respectively. | • | and 3f?(-) denotes the absolute value 
and the real part of a complex- valued number, respectively. ||a;||= V x^x and 1^ = denote the 
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vector 2-norm and vector direction, respectively. _L, || and denote perpendicularity, paralleUty and 
unparallelity, respectively. I denotes an identity matrix of an appropriate dimension. Rank(X) and 
Ti{X) denotes the rank and trace of the matrix X, respectively. X ^ 0, X ;^ 0, or X -< means 
X is a positive semidefinite matrix, a positive definite matrix, or a negative definite matrix. \i{X) and 
Ui{X) denote the i-th largest eigenvalue of X and its corresponding eigenvector, respectively. Xi{X, Y) 
and Ui{X, Y) denote the i-th largest generalized eigenvalue of X and Y its corresponding eigenvector, 
respectively. IIx = X(X^X)~^X^ denotes the orthogonal projection onto the column space of X, 
and = / — Ex denotes the orthogonal projection onto the orthogonal complement of the column 
space of X. 

II. System Model 

A. Signal Model 

Consider a two-user MIMO IC, where the transmitters and receivers are equipped with Nt > 2 and 
Nr > 2 antennas, respectively, and only one data stream is ttansmitted to each receiver. The data received 
by RXj is modeled as 

Ui = gif {HiiWiXi + HkiWkXk + ni), i, k e {1, 2} and k / i 

where Xi ~ CM{0, 1) is the transmitted symbol of TXj by the transmit beamformer Wi G C^^^^. At 
RXj, gij G C^"^^ is the receive beamformer, and rii G C^"^^ is the additive white Gaussian noise 
(AWGN) vector with zero-mean and covariance matrix afl. The matrices Hii,Hki € ([^n„xNt ^j-g jj^g 
channel-matrix of the direct link TXj i— > RXj and the cross-talk link TX^ H> RXj, respectively. Each 
transmitter has a power constraint that we, without loss of generality, set to 1 and define the set of feasible 
transmit beamformers as 

W = {«;GC^-^i:||w;f <1}. 

B. Rate with MMSE Receiver 

Assume that the interference from the other transmitter is treated as additive colored noise at each 
receiver. The rate of the link TXj i-> RXj is given by: 

Ri{wi,W2,gi) = log2 {l + SlNRi{wi,W2,gi)), (1) 
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where SINR4{wi,W2, Qi) = ^'^^^Ig^H^^Wkl'^ • linear transceiver design, it is well known that the 

MMSE filter is the optimal receiver for given transmit strategies. In this paper, the MMSE filter gij = 
{afl + HiiWiWi^ Hii^ + HkiWkWk^ Hki^y^ HiiWi is employed at RXj. Then, SINRi(it)i, iU2, fifj) 
becomes 



=Ai{wk 



SmRi(wi,W2) = wf Hfi{afl + HkiWkWk^Hki") ^HuWi. (2) 

In (O, the complex mathematical structure (inverse of the sum of matrices and product of matrices) 
causes a hard coupling of wi and W2 in the SINR in Q, which makes it difficult to analyze the SINR 
directly. To gain an insight into this coupling problem, we propose another form of the SINR expression. 

Proposition 1 For the two-user single-beam MIMO IC, the SINR in (jS]) can be reformulated as 

SmMwuW2) = sm\9H,) • "^"f"' + cos\9h,) • J'^""^^"' (3) 

af + \\HkiWk\\ 



where cos{6H,i) = \HiiWi ■ HkiWk\ and 9H,i & [0,tt/2]. □ 

Proof: See Appendix A. ■ 
Note that the SINRjfwi , W2) can be considered as a combination of ^^'^"^'^^ and H-^"'^-ll ^^jj-j, 
the weights sm'^{9H,i) and cos'^{9H,i)- That is, SINRj(i(;i, 102) depends not only on the power of the 
desired signal and the power of the interference HiJ^jiUfcH , but also on the Hermitian angle 



9H,i between the directions HuWi and Hj^iWf^. The SINR is coupled in a difficult way because of the 
existence of 9H,i- This is why it is more difficult to analyze the SINR of a MIMO IC than that of a MISO 
IC. How to decouple this hard coupling problem to simplify the SINR analysis is a main challenge in 
this paper. 

III. Pareto Boundary and Computation of Some Key Points 
A. Pareto Boundary 

The rate region is defined as a set of the achievable rate pairs with all the feasible beamformers 



■7^= [J iRl{wi,W2),R2{wi,W2)]. 
/Ill, ^ ' 



'Wi,'W2£W 

The outermost boundary of TZ is called Pareto boundary. The Pareto boundary can be divided by two 
turning points into the weak Pareto boundary (including the vertical part and the horizontal part) and the 
strict Pareto boundary (including the upper-right part and two turning points). 
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The Pareto boundary can be denoted by a set IZ* = U(^i' ^2) where {R\, R^) is any Pareto-optimal 
operating point. For {R\, R2) on the weak Pareto boundary, one rate is always fixed as the maximum rate. 
For (ii*,i?2) ori the strict Pareto boundary, it is impossible to improve one rate without simultaneously 
decreasing the other, and it can be reached unless both TXi and TX2 spend the maximum transmit 
power. More precisely, we have the following proposition. 

Proposition 2 For the two-user single-beam MIMO IC, all the operating points on the strict Pareto 
boundary can be achieved only when both the transmitters spend the full power, i.e., \\wi\\'^= ||w2p= l-D 

Proof: See Appendix B. ■ 
Although a similar result for the two-user MISO IC has been proved in |[T1, it cannot be extended to the 
MIMO case directly. As we will see later, this characterization of the strict Pareto-optimal beamformers 
is important for computation of the strict Pareto boundary. Here, we define a set of all the beamformers 
with full transmit power as 

W^P = {we C^-^^ : \\w\\^ = 1}. 

Note that all the strict Pareto-optimal transmit beamformers should be in the set Wj^-p. 

Due to the Pareto-efficiency, it is desired to find Pareto-optimal strategies to achieve the Pareto 
boundary. 

B. Computation of Some Key Points 

In this part, we compute exactly the Pareto-optimal points (the whole weak Pareto boundary and two 
turning points on the strict Pareto boundary) and some of the ZF operating points. 

1) Single-User Points SUl{Ri,0) and SU2{0, R2): A single-user point can be easily achieved when 
only TXj works and simultaneously operates "egoistically" to maximize its own rate. The maximum 
achievable rate Ri of the link TXj RX, and its associated "egoistic" strategy wf^° are 



i?i = Iog2(l + ^^^^^), (4a) 



'1 , ^ijHf^Hj, 

wf'° = ui{HfiHu), i = l,2. (4b) 

2) Turning Points Tl{Ri,R2) and T2{R^, R2): A turning point can be achieved when one transmitter 
employs an "altruistic" strategy to create no interference to the other receiver and simultaneously to 
maximize its own rate and the other transmitter operates "egoistically". For Tl{Ri,R2), we easily find 
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()h,i = 7r/2 in Q can result no interference in the cross-talk link TX2 1— ?• RXi. How to find the 
"altruistic" strategy ^2^'* is shown in the following proposition. 

Proposition 3 T'l(i?i,^2) ^'^'^ achieved by tu^'*), where Ri and w^^° are in (0) and 

R, = log, (1 + wf'^'A^iwfnwf') , (5a) 



= H.^f-^i n^£H..<-) ' (5b) 

Proof: See Appendix C. ■ 
Similarly, T2{R^,R2) with (wf * , W2^°) can be easily obtained by interchanging the indices. 
3) Weak Pareto-Optimal Points {R*, R,): For the weak Pareto boundary, both the horizontal part and 

the vertical part start end with a single user point and a turning point. Therefore, an arbitrary operating 

point {R\, R2) on the weak Pareto boundary can be computed as 

R^ = l-Ri and Rl = Rk, Vi, k£{l,2}, k^i 

where i = I and i = 2 correspond to the horizontal part and the vertical part, respectively. The scalar 7 
satisfies 7 G [0, 1). The operating point [R^jR,] becomes a single-user point or a turning point when 
7 = or 7 = 1, respectively. The associated weak Pareto-optimal transmit strategies are 

w* = ^-wf' and wl = w^^°. 



I 



Note that it is not necessary for both transmitters to spend the maximum transmit power simultaneously 
to achieve the weak Pareto boundary, which is different from the strict Pareto boundary (Proposition 2). 

4) Zero-Forcing (ZF) Points ZF{Rf^,R^^): In ©, we find 6h,i = 9h,2 = vr/2 results no 
interference in the cross-talk links TX2 1— )• RXi and TXi 1— )• RX2 simultaneously. The ZF conditions 
are 



> TT 

Gh,i = 7r/2 Hiiwi ■ H21W2 = 44> HnWi ± H21W2 <^ W2 1- H21H11W1 

>^ ^ N 

0H,2 = 7r/2 <^ H22W2 ■ H12W1 = 44- H22W2 -L H12W1 <^W2l- H22H12W1 
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from which and under a sufficient (and necessary only when Nt = Nj{ = 2) condition^, i.e., H21H11W1 \ 
H22H12W1, we obtain some of the ZF transmit strategies 

Nt-I 

wf^ = UiiHgHi2,H^iHu) and w^^ = J] CiU,{Ujj„j,^^^^,), 

1=1 

where {ci}fj^^^ are complex-valued numbers and satisfy Xli^i"^ kiP = 1- 
The ZF{Rf^,Rl^) can be achieved by (wf^,itjf^) as 

Rriwf^,wl^) = log, (1 + \\^^^^py , = 1,2 

Although the ZF points are not on the Pareto boundary, it is interesting to study ZF strategies if there 
exsits an additional requirement (like interference temperature or secrecy constraints) that each receiver 
is not allowed to receive the interference from other transmitters. 

IV. Computation of the Strict Pareto Boundary 

In this section, we propose an algorithm (called lAA) to compute the strict Pareto boundary. Since the 
rate region of the two-user single-beam MIMO IC is always a normal region according to Proposition 2, 
any operating point on the strict Pareto boundary is uniquely determined when one rate is fixed and the 
other rate is maximized. This motivates us to propose an optimization problem as follows 

max SINRi(iDi, lu,) 

(PO) < ■^1''^2GW^t> 

s.t. SINR2(wi,'u;2) = SINR2. 

where SINR^ e (2& - 1,2^^ - 1) is a SINR constraint, and w\^W2 should be in Wj-p according 
to Proposition 2. Then, an strict Pareto-optimal point (i?^,i?2) = (log2(l + SINRi(i(;^, lUg)), log2(l + 
SINR2(wf , 102))) can be achieved by the optimal solution (w\^w*2) to (PO). 

However, it is difficult to solve this non-convex optimization problem directly due to the hard coupling 
of tui and W2 in the SINR expression. To resolve this coupling problem, we convert (PO) into two single- 
beamformer optimization problems by fixing w\ and W2, respectively. Since both single-beamformer 
optimization problems are still non-convex, how to find the optimal solutions to these optimization 
problems will be studied in Section IV- A and Section IV-B. The Section IV-C will describe the algorithm 
in detail. 

'This condition is the same as that in 1151 , while we derive it from a different perspective (Hermitian angle in JS]!)- 
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A. Optimization of TXi 

For a given feasible W2 (how to find a feasible W2 will be shown in Proposition 4), the problem (PO) 
becomes a single-beamformer optimization problem with respect to wi. Its constraint is 

Hg{all + Huwiwi^ Hi2"y^ H22W2 = SINK* 

^W^HgH22W2 - l<^f2Jfl2t.lP ^ .gj^j,.^ 

C7| + ||iil2Wl|r 
(b)'Wi^H^2H22'W2wgHgHi2Wi HttHtt 2 



H^2H22W2 - cT^SINR^, (8a) 



J^w^C{w2)wi = and wgHgH22W2 > a^SINR^. (8b) 

The transformation (a) is based on the matrix inverse lemma. The transformation (6) is due to ||ij;i|p= 1. 
In the transformation (c), the nonnegative left-hand side of (H^) demands W2 H22H22W2 > (T2SINR2, 
and C is a Hermitian matrix defined as 

C{W2) = HgH22W2wgHgH^2 - {wgHgH22W2 - ajSlNRl) ■ {all + H^H 12)- (9) 
Then, wi can be optimized by 



(PI) 



max Wi Ai{w2)wi 



(10) 



S.t. Wi^ C{w2)wi = 

where C{w2) and Ai{w2) are Hermitian matrices. The problem (PI) is a homogeneous QCQP formu- 
lation, where both the objective function and constraints are quadratic without linear terms. It is difficult 
to solve this non-convex problem. 

Note that Xwi = Tt{XWi) for any matrix X, where Wi = wiw^ is a rank-one Hermitian 
positive semidefinite matrix. Thus, (PI) can be transformed to a convex problem by the SDR methods. 

' max Tt(Ai(w2)Wi) 

WitO V \ ' / 

(P2)<j s.t. Tt(c{w2)Wi) =0. (11) 
Tr(l^i) = 1. 

Observe that this convex problem is solvable (see the Appendix D). It is easy to obtain the optimal 
solution Wi by a convex optimization toolbox, e.g., SeDuMi ll20l . 
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Therefore, we can construct a tight optimal rank-one solution to (PI) by 

ui{Wl), if Rank(Vr^) = 1, 

wi = { (12) 
nV{W*^,C{w2),Ai{w2),l), if Rank(W^) > 2 and A^T > 3. 



where TZV denotes the function of matrix rank-one decomposition |[T9| . 
B. Optimization of TX2 

For a given feasible wi, (PO) becomes another single-beamformer optimization problem with respect 
to W2. Maximization of its objective function is 

max wfHfJajl + H2iW2W^H^A''^HiiWi 

_ I ( HrjHrr H^^H 2lW2\^ 

' — ^ max ■ [Wi HiiHiiWi 



af \ a\l + W2^ H21H21W2 

\W^H^^H21W2\'^ 



(13) 



mm 



lOsSW^p all + W2^ H21H21W2 
W2^Ci{wi)w2 



mm 



where Ci{wi) = HgHuWiwf H^^H2i and C2 = (rfl + HgH2i are Hermitian matrices. 
Then, W2 can be optimized by 

W2^Ci{wi)w2 



mm 



(P3) { '"^ew^p w^C2W2 (-24-) 
s.t. w^A2{wi)w2 = SINR^. 

Observe that (P3) is a fractional QCQP problem. It becomes the following SDR problem 

Tv{Ci{w^)W2) 



{P4){ 



mm 



W2^o Tr{C2W2) 

(15) 



s.t. Tr{A2{wi)W2) = SINR^ 
Tt{W2) = 1, 

which is still a non-convex problem. Fortunately, the fractional structure can be removed by a variation 
of the Chames-Cooper variable transformation |[2]1 . which has been used in |[22l . Define the transformed 
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variable Q = SW2 with s 




s.t. Tt{A2{wi)Q) = s ■ SINR^. 
(P5) I TV(C2Q) = 1 

TV(Q) = s 



(16) 




< s < . 

Ai(C2) AAr(C2) 



which is a standard SDR problem with two variables Q and s. We obtain the optimal solution W2 = 
to (P4) (see the Appendix E). 

Therefore, a tight optimal rank-one solution to (P3) is constructed by 




C. Algorithm 

Based on the above optimization, we propose an iterative algorithm to solve (PO) by optimizing wi 
and W2 alternately. We call this algorithm iterative alternating algorithm (lAA), which contains four steps: 
Step.l Set an initial feasible W2 for (PO). Step.2 Optimize (PO) with respect to wi by fixing W2. Step.3 
Optimize (PO) with respect to W2 by fixing wi. Step A Repeat Step.2 — 2> until this algorithm converges. 
As mentioned above, a feasible initiaUzation of W2 for (PO) is required, which will be discussed in the 
following. 

1 ) Initialization: A feasible initial W2 for optimization of TXi can be obtained as follows. 

Proposition 4 For a given SINR2 G (2—2 — 1,2^^ — 1), {wi,W2) is always a feasible solution pair to 
(PO) if Wi G yVjrp, W2 G Wjr with 



= \w2 eWTV-w^HgH22W2>a^SmR*2, Xi{C{w2)) ■ Xn{C{w2)) < o\ , (18) 



Note that Wjr is always nonempty because SINR2(iui, 1^2) = SINR2 with W2 G Wjr always 
corresponds to a straight line in the rate region. Therefore, the problem (PI), derived from (PO) with 



where C{w2) is in ((91). 



□ 



Proof: See Appendix F. 



February 27, 2012 



DRAFT 



13 



fixed W2 in Wj-, always has at least one solution wi in Wj-p by ([T2l ). 

Remark 1 worth to note that the performance of the lAA depends significantly on the initialization. 
Therefore, how to find a good initialization leading to a good/robust performance is studied as follows. 

In the two-user MISO IC, all possible strict Pareto-optimal beamformers can be characterized by linear 
combinations of ZF and MRT beamformers in jj^. A heuristical idea is to extend this parameterization 
for the MISO case to our MIMO case by balancing the "egoistic " strategy and "altruistic " strategy as 

W,{ii^,^^,2) = ^^,lwf'" + ^^,2wt^\ i = l,2, (19) 

where ^j^i and ^i^2 ^Jre complex-valued parameters satisfying i| + |^j^2| = 1- 

Unfortunately, the illustrations in Section V show that this characterization cannot achieve the whole 
strict Pareto boundary for the MIMO IC but still has a promising performance. To improve the efficiency 
of the initialization, we set an initial W2 by the following parameterization with real-valued parameters 

W2 = C-w^<^°+{l-C)-wf\ (20) 

where C = -p^ ~^ denotes the proportion of the ''egoistic" strategy. If W2 4. Wj-, we reset C, E + 

H2 — H2 /t2 — Hp 

— ^ -f , =^ — -2- J- until W2 G yVj^- If we still cannot find a feasible W2for all 
the V, we use a random W2 & Wjr directly. 

The later simulations in Section V imply this proposed initialization is a good candidate for the lAA. □ 

2 ) Convergence: We repeat the optimization of TXi and TX2 alternatively until a certain convergence 
criterion is satisfied, e.g., — < e with / and e as index of iteration and the accuracy, 
respecitvely. Since Ri is bounded by (0, Ri) and monotonically increases in each iteration, convergence 
of Ri is guaranteed. Finally, we obtain a convergent operating point iJg) with its associated transmit 
beamformer pair {wf\w2^). However, it remains unknown how close a convergent point (i?['\i?2) is 
to the strict Pareto-optimal point {R\, R2) (global optimality). Numerical simulations in Section V show 
that the lAA has a better performance to illustrate the strict Pareto boundary compared with the existing 
methods. 

3) Algorithm in Pseudo-Code: The outline of the proposed lAA is described in pseudo-code as 
Algorithm 1: 

4) Distributed Implementation: For the purpose of distributed implementation, we make two rea- 
sonable assumptions (similar to IIT2I). First, we assume that local CSI is available for each transmitter, 
namely, the transmitter TX^. knows the local channel matrices Hj^i to the receivers RXj,i = 1,2. The 
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Algorithm 1 Iterative Alternating Algorithm (lAA) 

Input: wfK an arbitrary R* € (^2>^2), and SINR^ = - 1. 

Output: A convergent operating point (i^^'^i^g) and its associated strategy {wf\w2^). 

begin 

Initialization: Set a feasible w'^^ G W^, / = 0. 
wliile some convergence criterion is not satisfied do 

Obtain an optimal solution VFi by solving the SDR problem (P2). 
if Rank(Vri) = 1 then 
I perform eigen-decomposition w^^ = ui(Wi); 
else 

|_ perform matrix rank-one decomposition wf^ = TZVfWi, Ai{w2~^^),C{w2~^^), l). 

Obtain an optimal solution (Q, s) by solving the SDR problem (P5). Set W2 = ^• 
if Rank(W2) = 1 then 
I perform eigen-decomposition W2'' = ui{W2)', 
else 

|_ perform matrix rank-one decomposition W2^ = Tl'D(W2,Ci{wf^),A2{wf^),C2,l)- 
\_ Compute the rate of hnk TXi ^ RXi: i??^ = logg (l + wf^'" Ai{w^2^)wf^). 



second assumption is that there exists an additional channel between the TXi and TX2 for information 
exchange. Under these two assumptions, the proposed lAA can be implemented in a distributed fashion. 

More precisely, for the optimization of TXi, TXi estimates the local channel matrices Hu and if 12 
and also needs the updated information (i.e., two vectors H21W2 and H22W2 and two scalars crl and 
SINR^) from TX2 via the additional channel between the transmitters. For the optimization of TX2, 
TX2 also needs the local channel matrices H22,H2i and the feedback (i.e., two vectors HuWi and 
H12W1 and a scalar af) from TXi. Furthermore, it is easy to stop the iteration when the termination 
criterion of the lAA is satisfied at TXi. 

Although the distributed implementation of the lAA needs information exchange, the exchange is only 
required between the transmitters (e.g., BSs for the downlink). This reduces the computation requirements 
at the receivers (e.g., MSs for the downlink). 

V. Illustrations 

To illustrate the achievable rate region of the proposed algorithm, we consider a two-user Gaussian 
MIMO IC, where Nt = 3 and Nr = 2. The transmit power budget is set to 1 for the two users, and 



February 27, 2012 



DRAFT 



15 



noise power af = a2 = 10 ") where SNR=10dB. The channels are 



H 



11 



Hi2 — 



H 



21 



-0.3034 + 1.9096i -0.3790 + 0.4201i 0.0357 + 0.7337i 

-0.6358 - 0.8030? -0.7881 - 0.1273f 0.7534 + 0.8348i^ 

-0.6758 + 0.1040? -0.5949 - 0.0344i 0.4311 + 0.9658i' 

-2.1621 + 0.5404Z -0.0037 + 0.6627? 0.8611 + 1.2318?^ 

0.3999 + 0.1567i 0.3798 - 0.5619i -0.1005 + 0.2836i^ 

-0.5494 - 0.4648? 1.1971 - 0.5297i -0.7271 + 0.2114f 



H 



22 



-0.0308 - 0.1133i 0.0433 - 0.3313i 0.3047- 1.2157? 
-1.4947 - 1.8676? -0.9430 + 0.5704? -1.3328 -M.4638? 



Fig. 1 shows the performance of the lAA by different initiaUzations. For a given R2 = R2 + {R2 — 
R2) = 6.2898 (far away from the turning points), Fig. la shows that the achieved Ri with initialization 
by ( [20I ) is greater than that with random initializations except for 8 out of 200 random initializations. For 



a given R2 = R2 + ■ {R2 — R2) = 5.6398 (close to a turning point), Fig. lb shows that the achieved 
Ri with initialization by (l20l) always outperforms the case with 200 random initializations. Therefore, 
numerical simulations show that the lAA with initialization by ( [20l ) has a good performance. 

Fig. 2 illustrates the achievable boundary by the proposed lAA compared with the existing methods. 
The term "IAA_10rand" denotes the best result obtained by running the lAA with 9 randomly chosen 
initiaUzations and 1 initiaUzation by (|20l). The SINR targets are SINR2 = 2^2s where = R2 + 
^ • (R^-Rg), n = 1,2, ...,49. The term "WMMSE_1 Grand" is defined with 10 randomly chosen 
initializations. The weighted sum rate used in "WMMSE_1 Grand" is expressed us, w ■ Ri + {1 — w) ■ R2 
with the weights ws, in [0.05 : 0.05 : 0.95]. The curve denoted by "RandBeam_lGmir' means the 
outermost boundary of the rate region achieved by 10 million random normalized transmit beamformers 
(Each receiver is fixed as the MMSE receiver). The term "SimpleReceiver" is the outermost boundary 
of the region achieved by ||5l where each receiver is fixed as the largest left singular vector of each 
corresponding direct channel matrix. The curve denoted by "Eq.([T9l)" is the outermost boundary of the 
achieved region by ( fT9l ) with complex-valued parameters. The boundary of "Eq.(l20l)" illustrates (l20l ) with 
Q = n = 1, N + 1 where N = 100. 
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It can be observed that the proposed lAA and the WMMSE yield a similar performance at convex parts 
of boundary and outperform the others under the same accuracy for convergence. However, since the 
WMMSE in 191 is a weighted sum-rate maximization method, it cannot achieve those operating points of 
the non-convex boundary. This might be why there is no operating points achieved between the points Pi 
and P2 by in Fig. 2. In addition, the performance of the "IAA_10rand" cannot be affected significantly 
by the number of random initializations from Fig. 1. 

To further evaluate the performance of "WMMSE_1 Grand" and "IAA_1 Grand" on illustrating the Pareto 
boundary, another simulation is done and shown in Fig. 3. Even with fine weights ws in [0.05 : 0.005 : 
0.95], we find there exists a large jump between the points P3 and P4 by "WMMSE_1 Grand" so that the 
rate region cannot be illustrated effectively. For "IAA_lGrand", although there is no guarantee to say that 
its corresponding curve is the strict Pareto boundary, it can serve as a good lower bound for the strict 
Pareto boundary. 

VI. Conclusions 

In this paper, we study an open topic on the Pareto boundary of the rate region for the two-user 
single-beam MIMO IC. The turning points and the weak Pareto boundary can be computed exactly. For 
computation of a strict Pareto-optimal point, our proposed lAA converges to an suboptimal operating 
point. The merit of this algorithm is that it can achieve any operating point with a given rate constraint. 
Thus, it can be appUed to the optimization problem with a rate constraint (e.g. rate requirement of primary 
link in underlay cognitive radio environments or for private messages), and it can illustrate the rate region 
effectively compared with the existing methods. 
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Appendix 

A. Proof of Proposition 1 

Proof: Given beamformers wi and W2, according to the matrix inversion lemma ||23l . Q can be 



rewritten as: 



SlNRi{wi,W2) = (Huwi) { —I 2/ 2 ■ "^^ —]Huwi 



al al{al + \\H2iW2\\'^] 
\HiiWi\\^ \{HiiWi)^H2iW2\ 



o-f a{{a{ + \\H2iW2\ 



\HiiWi\\^ Hiiw\^ ■ H2iwl ^ • ||i3"21''A'2| 



a( V af + \\H2iW2\\^ / 

T^H 77 ^|2\ ||i3"iiit)i||^ |— ,„> |2 ||i3"iiiui|'^ 



1 - i3"lltOl • i3"2lt«2 • 9 h i3"llWl • i3"21'«A'2 



0"l erf + ||i3"2llf2| 

(21) 

For two complex vectors a and b, the cosine of the complex-valued angle between a and b is defined 
as Il24l: 

a^b 



cos(^i 



||a||-||&|| 

where 003(6*0) = /^e-^'^ with ^ = | cos{0c)\ < 1 and — vr < < ^ is called pseudo angle between a 
and 6. 

The Hermitian angle between a and b is defined as 

\a^b\ 

cos(0h) = I cos(%)| = ' ,, < < tt/2. 

||a||-||o|| 

y ^2 ^ii2 ^ii2 

It implies \HiiWi^ ■ i3"2i''A'2| = cos'^{9h,i) because of = ||ii"2iii'2|| = 1- Thus, (|2TI ) 

becomes 

1 1 2 1 1 2 

SINRi (^1,1.2) = sm\eH,i) ■ "^^T^" +cos^(gH.i) • ^ n^'J^'" ,,2 , 

O"! Cjf + ||i?2l1t>2|| 



where 9h,i G [0, vr/2] denotes the Hermitian angle between the desired signal direction HnWi and the 



interference direction H21W2 at RXi. Obviously, when HnWi \\ H21W2 (or HnWi ± H21W2), we 
have 9h,i = (or 9h,i = 7r/2). ■ 
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B. Proof of Proposition 2 

Proof: Our idea is to show that it is impossible for an strict Pareto-optimal point achieved by the 
transmit beamformers with less than full power. Therefore, the proof works by contradiction. 

Assume that an operating point W2), R2{wi,W2)) on the strict Pareto boundary is achieved by 

||wip< 1 and ||i02p< 1- We consider whether there exists an outer point (^Ri{wi,W2), R2{'Wi,W2)) 
achieved by \\wi |p< ||'u;ip< 1 and ||i02|p< 1. If an outer point exists, e.g., Ri{wi, W2) > Ri{wi, W2) 
and R2{wi,W2) = R2{wi,W2), we can improve the Ri{wi,W2) only by more transmit power while 
keeping R2{wi,W2) unchanged. Thus, the existence of an outer operating point contradicts the assump- 
tation that (^Ri{wi,W2),R2{wi,W2)) is a strict Pareto-optimal point. 

That is, we need to show the existence of a nonzero perturbation vector dp satisfying: 



{wi + Sp)^ Ai{w2){wi + Sp) > wfAi{w2)wi (23a) 

W2A2{wi + Sp)w2 = W2 A2{wi)w2 (23b) 

\\wi + 5p|p> \\wi\\'^ (23c) 

\\wi + Spf<l. (23d) 

An arbitrary nonzero Sp can be expressed as 

Sp = \\dp\\-e^^' -Sl. (24) 



It means that we should find ||^p||, (ps and Sp to satisfy all the conditions in ( |23] ) simultaneously. The 
proof of the existence of (ps is similar to that for the two-user MISO IC in [1]. However, it is more 
difficult to find a Sp for the MIMO IC because the cross-talk channel matrix (rather than a vector in the 
MISO IC case) does not always have a null space for S^. We give the proof in detail as follows. 
1. Existence of S^ 

By the matrix inverse lemma ll23l . the condition in (l23b) is equivalent to 

\{wi + Sp)^H^2H22W2\^ ^ \W^H^2H22W2\^ ^^5) 

a^ + \\Hi2iwi + Sp)f a^ + \\Hi2Wif ' 
It is difficult to solve Sp directly. In fact, we only need to prove the existence of Sp satisfying (1251 ). 
Case 1. Nn < Nt or Rank(ii'i2) < Nt < Nr: 
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We always have Hi2Sp = if 

Rank(n-L ^ ) 

Sp= ^ aiUi{UjjTj, (26) 

Rank(n^y ) 

where ai,i = 1, Rank(n^T ) are complex-valued numbers and J2i=i = 1- Then, ( |25l ) 

always holds because any 6p in the null space of the cross-talk channel H12 does not cause extra 
interference to RX2. 

Case 2. Rank(i3"i2) = Nt < Nr: 

It is impossible to nullify the perturbation directly as Case 1. Define vi = H12W1, vs = Hi2^p and 
V2 = H22W2- Then ( |25] | becomes 

9ll ||2 9iiii2' v^'/ 

+ ||f 1 + f^ll O'l + ll'l'lll 

Assume that is a combination of two orthogonal vectors 

"^5 = n^^i + - ^ • '^v^Vl), (28) 

Note that H^jt'i = V2 ■ e~^'^^ where (pi = aTg{v^V2). Now, it remains to find whether there is a vs in 
the plane spanned by H^^i'i and fUt^i satisfying (|27] |. 
Substituting ([28]) into ([27]) yields 

l^f + ||^'<5|| • yi - T? • e^*?^^ • \\V2\\\^ _ |t7ft>2p 

^ 2 2 II II 2' v^"/ 

+ + • {Vv ■ ^i^i + ■ e-^^^ ■ ^) || ^2 + ll^ill 

Define the right-hand side and the left-hand side of ( [29l ) as Rside and Lsideiv)' respectively. It is still 
hard to get a closed-form solution of rj by solving ( |29l ) directly. Observe that the denominator of Lgide (v) 
is always positive for t] G [0, 1], and Lgideif]) ^ function of t] is continuous over the interval [0, 1]. 
Therefore, if {Rside - Lside{l)) {Rside - Lside{0)) < 0, there must exist a vs{'n) with some t] G [0, 1] 
satisfying ( [291 ). 

When 7] = I, the term Lside{f]) becomes 

T n\ - \viV2\'^ _ \vfv2\^ 

J^side{^) — — 



£72 + + llf^ll • n^^f ill 2 II i an II V2 il 

Observe that WvsW^ + 2||d5|| • ||n^ i'i||> 0. Thus, we have Lside{'^) < Rside- 
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When 7] = 0, the term Lsideirj) becomes 



\v^V2 + \\v5\\-e^'^- ■\\V2\\\^ ^ \v'^V2? + \\v2fi\\v5\? + 2\\V5\\-K^2\) 
al + 11^1 + ll^^ll • e-^^^ ■ ^2f ai + ||^iP+(||^5P + 2||t,5||-K^^|) ' 



where HiJ^p + 2 • \\v^\\ ■ \v^V2\ > 0. If Lside{0) > Rside, we need ||^'2|P> Rside- Furthermore, R^ide is 
bounded by 

_ |t'ff2p , ll^'l|P-|l^2|P ||f2|P II ||2 

^^'^^ ~ .2^11,, ||2 ^ ^2^11,, ||2 - " < ll'"2ll • 

ii^iir 



Thus, we have Lside(O) > inside- 

Due to Lside{l) < Rside < Lside{0), there exists at least one rjo G (0, 1) satisfying Lsideim) = Rside- 
In this case, H12 has an inverse/Moore-Penrose pseudo-inverse matrix h\2. Then, we have 



, (32) 
= \\8p\W'i'^ -5^, 



1 1 (5 II 

where ||t'5||= -n — -7 > s ^ -TT■ depends on but has no requirement for \\8p\\. Therefore, 

||-H"l2(,v^-n^2'"i+Vi-^o-n^2'"i ) 



vs with any \\5p\\ and ^ = e ■^'^'^ • i?|2(\/% " ^V2'"^ + \/l — r/o • n^j^i^i) satisfies (l27l) . 
Therefore, any 6p = ||(5p||-e-''^'' • 5p with 

Rank(n^ ^ ) Rank(n^ r ) 

^12 ^12 

^ aitti(n^T ), ^ |aip = 1, when iV/f^ < Nt or Rank(iJi2) < Nt < Nr 



dp 



i=l i=l 



^ e-^'^^ • J/la (\/% • n^^i + • n„,7;i) , when Rank(l/i2) = A^r < iV/j, 



(33) 



satisfies the condition (l23b). 
2. Existence of (ps 
Substituting dMll into (l23b) yields 



(Wl + 5p)^Ai(w2)(l(^l + 5p) > W^Ai{w2)Wi 

\\Spfdl''Ai{w2)t + 2||5pP(tz;f Ai(t«2)? • e^<^^) > 
^^dJ^Ai(i(;2)^ + \w^Ai{w2)sl\ ■ cos(</.5 + </>2) > 

^-^-N^ ll^pll sl"Ai{w2)sl 
cos((/>5 + 02) > • - — (34) 

^ Ai(w2)<5p 



February 27, 2012 



DRAFT 



21 



where (/)2 = arg{wi Ai{w2)Sp). 

At the same time, substituting (|24l) into (l23b) yields 

llifi + ^pll^> ll'ifil 



^ ||<5pf+2||dp||5?(t(;f5j-e^'<^^) > 
|i(;f ^1 cos(</>5 + </>3) > --^^^ 

II ^ II 

^ cos(05 + 03) > (35) 

2|wf5p| 

where 03 = arg(i(;f 5p). 

Define 05 + 02 G [6*1, 6*2] and 05 + 03 G [6*3, 6*4]- Since both the right-hand side of and (1351 ) 
are negative, the range [^1,^2] and [^3,^4] are strictly wider than vr. In addition, the intersection of two 
angular ranges wider than n is nonempty. Then, any Sp in (l24l) with 0^ G [^1,^2] H [^3,^4] satisfies the 
conditions (l34l) and (l35l) simultaneously. 

3. Existence of \\Sp\\ 

The condition (l23ll) is equivalent to \\wi + Sp\\< 1. In addition, due to + (5p||< ||wi|| + ||<5p||, any 
dp in dull with 

||^p||< 1 - Ikill (36) 

satisfies the condition (l23h). 

Above all, the existence of Sp, 4>s and \\Sp\\ has been proved. That is, there always exists some 
5p = \\5p\\-e^'^'^ ■ Sp satisfying all the conditions in (l23l) . Then, Ri{wi,W2) can still be improved untill 
||i(;i|p= 1, while Ri{wi,W2) remains unchanged simultaneously. This contradicts the assumption that 
[Ri{wi,W2), R2{wi,W2)) is on the strict Pareto boundary. Therefore, Proposition 2 holds. ■ 

C. Proof of Proposition 3 

Proof: For the turning point Tl{Ri, R2), to achieve the maximum rate of link TXi 1— RXi, i.e., 
Ri in (|4^), (11)1,102) should satisfy the following conditions 

= wf'° = «i(i/fi/fn), (37a) 
9H,i = 7r/2^W2±H^iHuWi. (37b) 
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This means W2 should be in the null space of Ht^iHuw^^" to create no interference to RXi. All 
If 2 G Wj-p satisfying (l37l) form a set of the ZF strategies, defined as WzT- Then, any W2 G W^j:- can 
be expressed as 



^2 = n^„^^^^.,„^2, (38) 

where V2 G C^^^^ and # /ffi/Zniff^". 

To achieve (i?i,i?2)> we need to find i^g^* which maximizes SINR2(iyf 1^2) simultaneously. Here, 
we define the optimal "altruistic" strategy lUg^'* as 

lu^'* = arg max W2 A2{wf^°)w2 

^lPt,f = arg max ^ j^fii^..^- ^ ^ m.H^.^f^^ 



^2 nH-H,,^f-^2(ti^f ^")n^„ V2 

= arg max \^ — j 

= «i(^i'n^-i^.<-)- (39) 

Substituting ( [39] ) into ( [38] ). we obtain the optimal "altruistic" strategy 



Therefore, Tl(Ri , R2) can be computed as ( logs (} + ^'^•^jj^"^ ) , loga (l + wf*'^ A2{w^^")wf^ 
by (tff^^tff*). ■ 

D. Solving the SDR problem (P2) 

Before we solve the problem (P2), we need first to consider its solvability. 

Lemma 1 The problem (P2) is solvable, i.e., its respective finite optimal solutions exist. □ 

Proof: The given W2 is assumed feasible for (PO). It means that there exists at least one solution 
wi € yVj^-p to (PI) by fixing W2. Thus, Wi = wiw^ should be a solution to (P2) and its objective 
function is still bounded by (0, 2^^ — 1). Based on Weierstrass' Theorem, the problem (P2) is solvable, 
i.e., its respective finite optimal solutions exist. ■ 
We can obtain the optimal solution W\ to the problem (P2) by a convex optimization toolbox. 
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However, the rank of an optimal solution to (P2) is usually more than one because we have discarded 
the rank constraint Rank(Vl^i) = 1. If W\ is rank-one, then we can easily get the optimal rank-one 
solution to (P2) by eigen-decomposition, i.e., wi = ui{Wi). If the rank of Wi is greater than one, we 
can perform a matrix rank-one decomposition W* = J2^=i^^^^'^ x*x*^ such that 

X, C{w, )x,- j^^^j^^^.) and x, - ^^^^^^.^y Vr - 1, Kank(M/J. 

Since (P2) has two constraints, there is no gap between the optimal value of the SDR problem (P2) of 
X* and that of the QCQP problem (PI) of wi [25|. The theorem and algorithm of the matrix rank-one 
decomposition are based on the results in |[T8]| . |[T9]| . Here, we quote Theorem 2.2 in |[T9]| in the following 
lemma. 

Lemma 2 Let Ai, i = 1,2,3 be N x N Hermitian matrices, and X be a nonzero N x N Hermitian 
positive semidefinite matrix where N > 3. Then, 

Tjf Rank(X) > 3, one can find in polynomial-time a nonzero vector y G Range(X) such that 

Tr{Aiyy^) = Tr(AiX), i = 1, 2, 3. 

^^^^ ^ - Rank(X) ^^^ ^ ^""^ Rank(X - Rank(X) ^^^) ^ R.ank(X) - 1. 

Tjf Rank(X) = 2, then for any z ^ Range(X) there exits one y G in the linear subspace spanned 
by zL) Range(JC), such that 

TriAiyy") = Tr(AiX), i = 1, 2, 3. 
with X + zz^- ^^yy"" h and Rank(X + zz^ - ^^^yy"") < 2. □ 

The only condition for the applicability of Lemma 2 to W* is Nt > 3. It means TXi should has 
Nt > 3 antennas. Then, we can obtain an optimal rank-one solution wi to (PI) via solving (P2). 

E. Solving the SDR problem (P4) 

We show the SDR problem (P4) and the fractional SDR problem (P5) are solvable and equivalent. 

Lemma 3 Both the fractional SDR problem (P4) and the SDR problem (P5) are solvable. □ 

Proof: For the fractional SDR problem (P4), its feasible region 

n = {W2 h : Tt{W2) = 1, Tj:{A2{wi)W2) = SINR^} 
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is a nonempty and compact set because both (P4) and (P2) are transformed from (PO). If Wi is a 
solution to (P2) for a given W2, W2 should be a solution (regardless whether it is optimal) to (P4) for 
a given Wi. 

Observe that both the denominator and numerator of the objective function of (P4) are continuous 
over the region Q. The denominator is positive and satisfies 

aj + XNAHgH2i) = XnAC2) < Tv{C2W2) < Xi{C2) = + Xi{HgH2i), (40) 

and the numerator is obviously nonnegative and satisfies 

< Tt{Ci{wi)W2) < Xi{HgHuH^iH2i). (41) 

This implies that the objective function of the fractional SDR problem (P4) is bounded over the region 
n by 

^ ^ Tt{Ci{wi)W2) ^ Xi{H^,HuH?iH2i) 



T^{C2W2) - al + XNAH^iH2i) ' 
Based on Weierstrass' Theorem, the problem (P4) always has an optimal solution. 

Assume that W% is an optimal solution to (P4), we know that s* = — 7—^ ^ and Q* = s* ■ 

are feasible for (P5). Also note that its objective function is also bounded by (|4T| |. Similarly, (P5) is 
solvable according to Weierstrass' Theorem. ■ 

Lemma 4 The problems (P4) and (P5) have the same value. Furthermore, if W\ solves (P4), then 
s* = —r^ ^ and = s* solves (P5); if Q" and s* solve (P5), then W% =% solves (P4).n 

Proof: Assume that W2 is an optimal solution to (P4), and v*p^^ and v*p^-^ are the optimal values 
of the objective function of (P4) and (P5), respectively. Thus, s = ^ ^'^'^ Q = s ■ W2 are 

feasible for (P5). The value of the objective function of (P5) at this feasible point is 

V(^P5)=Tr{Ciiwi)Q) 



Tv{Ci{wi){s-W'2)) 



Tr{Ci{wi)Wl) 
Tt{C2W*2) 



V 



{P4) 



On the other hand, suppose that Q* and s* are the optimal solutions to (P5). Since s* is always 
positive according to (|40l ). W2 = ^ is also feasible for (P4). Then, the value of the objective function 



February 27, 2012 



DRAFT 



25 



of (P4) at this feasible point is 

Tt{Ci{wi)W2) 



-'(P4) 



Tt{C2W2) 

Tr(Ci(t«i)2:) _ Tt{Ci{wi)Q* 



TT(Ciiw,)Q* 



\P5) 



Above all, we have v^-o.-, = 

{P4) (P5) 



Therefore, we can solve the fractional SDR problem (P4) via the solution to the SDR problem (P5). 
After obtaining the optimal solution (Q*, s*) to (P5), we get the optimal solution to (P4) by W2 = ^ 
based on Lemma 4. If W2 is of rank one, we get W2 by eigen-decomposition. Otherwise, we can 
construct a optimal rank-one solution W2W2 = W2 by a matrix rank-one decomposition in |[T9l . which 
also certifies (P5) is tight. Here, we quote the following lemma from the Theorem 2.3 in (T9\. 

Lemma 5 Let Ai, i = 1,2,3,4: be N x N Hermitian matrices, and X be a nonzero N x N Hermitian 
positive semidefinite matrix where N > 3. Furthermore, suppose that 

Tt{AiY), Tt{A2Y), Tr(A3l^), Tr(A4l^)) / (0, 0, 0, 0) (42) 

for any nonzero N x N Hermitian positive semidefinite matrix Y. Then, 

Tjf Rank(X) > 3, one can find in polynomial-time a nonzero vector y G Range(X) such that 

Ti{A,yy^) = TY{AiX), i = 1,2,3,4. 

7/'Rank(X) = 2, then for any z ^ Range(J^) there exists one y S in the linear subspace spanned 
by zL) Range (JC), such that 

Ti{Aiyy^) = Tr{AiX), i = l,2, 3, 4. □ 

It remains to show whether the Lemma 5 is applicable to the solution W2 to (P5). One condition 
Nt > 3 means that TX2 should be equipped with Nt > 3 antennas. The other condition in (l42l ) suffices 
to prove that there is (ai, 02, 03, 04) G M'^ such that aiCi{wf + a2^2(i'^i + «3C'2 + o-aI >- 
|[22ll . which can be satisfied by (ai, 02, as, a^) = (0, 0, 0, |a4|) with 04 7^ 0. 
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F. Proof of Proposition 4 

Proof: We need to find a feasible set Wjr such that there exists some solution G Wj-p to problem 
(PO) by fixing W2 G Wj-. 
In (Hh), the constraint of (PO) is equivalent to H22H22W2 > ct|SINR2 and 

w^C{w2)wi = 0. (43) 

To guarantee the existence of tui e Wj-p, a feasible W2 should be determined in a way such that ^p) 
is satisfied. 

Define C{w2) = Ylfji Xi{C{w2))ui{C{w2))uf {C{w2)). We analyze C{w2) for two cases. 

Case 1. When C{w2) is a full rank matrix, i.e., Xi{C{w2)) 7^ 0,Vi = 1,...,Nt, C{w2) has three 
possibilities: a positive definite matrix {C{w2) >~ 0), a negative definite matrix {C{w2) -< 0), and 
Xi{C{w2)) > and Xn^{C{w2)) < 0. 

If C{w2) is a positive/negative definite matrix, it is clear that there is not a nonzero vector W2 satisfying 
(|43] ). For C{w2) with Xi{C{w2)) > and AAr^(C(i(;2)) < 0, a sufficient solution to (|43]l is 



/ XNTiC(w2)) . / Xl(C(w2)) 

= V a,(c(^,))-a..(c(^.)) ^ • + V a,(c(^.))-a..(c(^.)) • --^(^(-^))- 



Case 2. When C(i(;2) is not a full rank matrix, i.e., Xi(C{w2)) = for some i G {1, Nt^, C{w2) 
has three probabilities: X\(C(w2)) = 0, Xj\[r^{C{w2)) = 0, and Ai(C(i(;2)) = for some i 7^ {1,Nt}- 
In this case, C(i(;2) always has null space for W2 to satisfy (1431 ). 

Above all, the sufficient and necessary condition of W2 satisfying (1431 ) is Ai (C(iU2)) • Xnj. [C{w2)) < 
0. Therefore, any feasible W2 for ([8J3) should be in the following set 

= jiua € W^v ■■ w^HgH22W2 > aiSINR^, Ai(C(k;2)) • Xn^{C{w2)) < 
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February 27, 2012 



DRAFT 



30 



<D 5.5 

CO 



CD 

c 
c 
CO 



tr" 4.5 













1 1 1 1 

- e - IAA_1 Grand 

- *■ - RandBeam_10mil 
< WMMSEJ Grand 
^ Turning point 

Y Turning point T^ 








Q 




• 

-te 


















V 






















1 










i 







2.5 3 3.5 4 4.5 5 5.5 6 

[bits/channel use] 



Fig. 3. Achievable boundary witli SNR=10dB and Nt ^3,Nr = 2 



February 27, 2012 



DRAFT 



